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Abstract 



type. Our method is based on Newtonian potentials through the framework set forth in the case of 
dimension two from complex analysis point view in our recent work. 

1 Introduction 



>• . In this paper we consider solvability problem for nonlinear partial differential systems in R n . 

In dimension two, the same study was carried out in [PI] using complex analysis with high 
j^. ■ order Cauchy-Riemann operators. Here in dimension greater than two, we will solely base on 
CO ■ Newtonian potential and its Holder estimate in R n . We obtain existence of non-radial classical 
\Q ', solutions of nonlinear systems of Poisson type under rather general conditions. The following 
^ are the main results of this paper. 

Theorem 1.1. Let a(x,p,q,r) = (ai(x,p, q, r), ajsf(x,p, q, r)) be of class Cf oc (0 < a < 1), 
where x G R n ,p G R N , q G lR n <g) R N , and r G Sym(n) <g> R N . Assume that 

G (°) = °. C 1 ) 

V r a(0) = 0, (2) 
V^a(0) = 0. (3) 

Then the following system: u(x) = (ui(x), ...,Un(x)) : {\x\ < R} — > R N , 

Au(x) = a(x,ti(s) ) Vu(i),V 2 n(x)) (4) 

has infinitely many solutions of C 2+a ({\x\ < R}) of vanishing order two at the origin for suffi- 
ciently small values of R. Moreover, all these solutions are not radially symmetric. 

Here we make comments on notations used above. First the notation V r means that deriva- 
tives are taken with respect to the variables of r, and Sym(n) denotes the set of n x n symmetric 
matrices. A solution u is said of vanishing order two if u(0) = 0, Vtt(0) = 0, but V 2 w(0) ^ 0. 
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The vanishing order two ensures that the solutions obtained are non-trivial; namely they are 
neither constant nor linear ones. More importantly, the vanishing order two ensures that solu- 
tions of such are not radial, and this shows that the solutions could not come from solving an 
ODE if the equation could be ever possible. 

If a is independent of r, then the following existence theorem can be regarded as a variant of 
ODE with initial values. We point out, though, that there is no uniqueness result as in ODE; 
in fact instead there are infinitely many solutions by construction. 

Theorem 1.2. Let a(x,p,q) = (ai(x,p, q), a^(x,p, q)) be of class C^ 01 (1 < k < oo,0 < 
a < 1), where x G R n ,p G R N , and qeR n ® R N . Then, for any given c G R N , c x G R n <g> R N , 
the following system: u(x) = (ui(x), Un(x)) : {\x\ < R} — > R N , 

Au(x) = a(x,u(x),Vu(x)) (5) 
«(0) = c (6) 
Vu(0) = ci (7) 

has infinitely many solutions of C k+2+a ({\x\ < R}) for sufficiently small values of R. In partic- 
ular, all solutions are not radially symmetric. 

If a is independent of x, that is, the system is so-called autonomous, then we can solve 
non-trivial semi-global solutions, i.e., solutions that are defined in any given ball < R}. 

Theorem 1.3. Let a(p,q,r) = (ai(p, q, r), a^(p, q, r)) be of class Cf oc f20 < a < 1), where 
p G R N , q G R n (8) R N , and r G Sym(n) ® R N . Assume 

a(0) = 0, (8) 
Va(0) = 0. (9) 

Then the following system: u(x) = (ui(x), ...,un(x)) : {\x\ < R} — > R N , 

Au(x) = a(u(i),Vu(x),V 2 u(x)) (10) 

has infinitely many solutions in {\x\ < R\x G M n } of C 2+a with vanishing order two at the origin 
for any given value of R. Consequently, all these solutions are not radially symmetric. 

We remark that the conditions (8) and (9) imply that the nonlinearity of a does not contain 
linear terms and that u = is a trivial solution of (10). These conditions are necessary due to 
a well-know result of Osserman [O]. In fact, let us solve the scalar equation 



2-u 



Au = e 

with the initial value u(0) = a, which exists for a small R according to Theorem 1.2. On other 
hand, by Osserman's theorem applied to this case, we have 

1 

R < 



oU(0) 



Letting a — > +oo, we see that R — > 0. Of course, a(p) = e p does not satisfies the condition 
(8). This example also shows in Theorem 1.1 and 1.2, R needs to be small in general. Another 
example, we may consider, is the eigenvalue equation 

Au = Aw, 
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which, of course, has no non-zero solutions for most values of A. 

A classical and well-known result of Gidas-Ni-Nirenberg [GNN] says that if u is a positive 
solution of the Poisson equation 

Au = f{u) 

in \x\ < R with / C 1 and 

U\{\ X \=R} = 0, 

then u is radial. However their theorem does not prove the existence of such a solution. As an 
application of our results, we prove the following corollary for existence of non-radial solutions. 

Corollary 1.4. Let f(t) be a function of class C k+a (R) (k > 1,0 < a < 1). Assume that 
f(0) = /'(0) = 0. Then the equation 

Au = f(u) 

has, for any given R, infinitely many solutions of class C k+2+a ({\x\ < R}) which are of vanishing 
order 2 at the origin and are neither radially symmetric nor positive. 

For example, much study has been conducted on seeking positive solutions of the equation 

i i "+ 2 
Au = c\u\ »- 2 

in R ra . Corollary 1.4, however, implies that there are indeed infinitely many non-radial and 
non-positive solutions for a ball of any radius. However it does not conclude the existence in 
the whole M ra , which is of course in general impossible due to [O]. 

We would like to point out that for solvability of linear partial differential equations, there 
is a well-known so-called Nirenberg- Treves conjecture. This conjecture was recently solved by 
Dencker [D], following previous important works in [L], [NT], [H], [BF], [LE]. Our consideration 
on nonlinear cases is very different from linear ones technically. 

In the paper [PI] dealing with dimension two, complex analysis allows us to prove similar 
results with power of Laplace. In higher dimension, we have to overcome more difficulty in order 
to carry out the method of this paper and that of [PI] for the general system 

A m u(x) = a(x, u, Vu, V"V). 

The complete work will be in a future paper [P2]. 

This paper is organized as follows. First we introduce a Banach space with vanishing order 
from which we are seeking possible solutions. Then, we study Newtonian potential as an operator 
on the Banach space, which seems to have been overlooked in the literature. Finally, we construct 
an operator map from which we will try to produce a fixed point on a closed subset of the Banach 
space. A large portion of the paper is on estimating Holder norm of involved functions in order 
to apply Fixed point theorem. 

2 Function spaces and their norms 

In this paper throughout, we let D denote the closed ball {x G M n | \x\ < R} and C its boundary 
{x G I" \x\ = R}. Unless otherwise stated, all functions considered will be real- valued and 
integrable, with domain D. We will consider some classes of functions. 
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2.1 Holder space 



C a (D) is the set of all functions / on D for which 

!/(*)- /(*')! 



H a \f] = SU P j" 



x, x' G D 



is finite. 

C k (D) is the set of all function / on D whose k th order partial derivatives exist and are 
continuous, k an integer, k > 0. C k+a (D) is the set of all functions / on D whose k th order 
partial derivatives exist and belong to C a (D). 

The symbol |/| or |/| D denotes sup xeD \f(x)\. For / G C a (D) we define the norm 

11/11 = l/l + ( 2 -R) Q #a[/]. 

The set of iV-tuples f — (f ± , f N ) of functions (vector functions or maps) of C a (D) is denoted 
by [C a (D)] N , and H a [f] is defined as the maximum of H a [fi](i = 1, .., N). In a similar fashion 
we define \/\a = sup x( z A \f(x)\ for functions and vector functions, and write |/| when the domain 
is understood. Finally, in this paper throughout, the norm of M. N is taken as \v\ = maxi<j<jv \vj\. 
The following lemma is well-known; see ([GT]). 

Lemma 2.1. The function || • • • || defined on C a (D) is a norm, with respect to which C a (D) is 
a Banach algebra: \\fg\\ < ||/||||^||. 



2.2 Function spaces with vanishing order at the origin 

Our idea of solving differential equations or systems of order m is to look for solutions that 
vanish up to m — 1 order at the origin; this way the norm estimate of the function space to be 
considered later is made possible using only mth order derivatives. This is rather different from 
classical norms used for higher order derivatives in partial differential equations. The results 
in this section have been given in [PI] in dimension two using all complex notations, and the 
method would work for any dimension. For convenience of readers, we produce all results in 
real variable notations in IR n . 

We denote for k > 1, Cq + °(D) the set of all functions in C k+a {D) whose derivatives vanish 
up to order k — 1 at the origin. Specifically 

C k+a (D) = {fe C k+a (D)\d?f(0) = O,\0\<h- 1}, 

where we have used /? = /?„) and \/3\ — j3i + ... + (3 n . Also we have 

One has the following obvious nesting 

C™ +a (D) C C™~ 1+a {D) C • • • C Cl +a {D) C C a {D). 
We now define a function || • • • ||^ on C k+a (D): 

||/||« = max{||^/||}. 

\B\=k 
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We point out that the function || ■ ■ ■ ||^ on C k+a (D) is not a norm since ||/||^ = if and only 
if / is a polynomial of degree at most k — 1. However it becomes norm when restricted to the 
subspace Cq +c *(D), which is to be proved below and is one of important facts used in this paper. 
First we obtain some useful estimates. 

Lemma 2.2. If f e C k+a (D), then, forx,x' e D, 

- E \ E d "f^ x ' ~ X ) P \ ^ hi E H ^f])\x' - x\ k+a . 
1=0 ' |/3H ' L W=k } 

Proof. Expanding at x, we have the formula 

fc-i , 

/(^-EiyE^)^'-^ 

Z=0 ' |/3|=i 

= l £ 1 "' r {^ft tx ' + ( i - t ^} dtdti --- dtk - 1 

= I f" \.. f 1 | J2 d fi f(tx' + (l-t)x)(x'-x) f) \dtdt 1 ---dt k - 1 . 
Jo Jo Jo I J 



|/3|=fc 

Hence, we have, by subtracting kth term, 

k 



/oo-E^E^x*'-*) 



Z! 

z=o 



1 r t 



'0 ./0 

Thus we have, 



{d^f(tx' + (1 - t)x) - d p f(x)}(x' - xf \dtdt! ■ ■ • dt fc _i. 



fc 

/! 

1=0 \/3\=l 



i/^-E^E^x*'-* 



<^E^[^]i x/ - x i fe+a 



fc! 

<8=* 



This completes the proof. 
Lemma 2.3. // / e C% +a {D), then 



< ^^ll/ll (fc) - 
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Proof. Let / G C% +a (D), then 



1 rtk- 



k-l 



r d k 



f(x) = I J ■ ■ I {-^f{tx)}dtdt 1 ---dt k - l 



o \dt k 

Jo Jo Jo l m=k J 

^ ••• jf & 3 f(tx)dtdt 1 ---dt k _ 1 j x 



/3 



Applying norm inequality, we obtain 

,/3 



(fc) 



< £^ii#7iiik\. 

|/9|=fc 
|/3|=fe 

where we have used = 3/2, which is easily verified. 
Lemma 2.4. If f e C k+a (D), then, for \/3\ < k, 

Proof. Let / G C fc+Q (D). If \f3\ < k, then d?f G C M/3|+Q? (D). By Lemma 2.3, we have 



An immediate corollary is the following: 
Lemma 2.5. If f E C% +a (D), then, for I < k, 

ll/ll (0 < j^R k - l \\f\\ ik) . 

In order to verify that Cq +(X (D) is a Banach space with norm || • • • ||^, we need the following 
simple lemma similar to ([NW], 7.1. a). 

Lemma 2.6. Let {fm\m=i ^ e a sequence in C a (D), with \\f m \\ < M; assume that {f m } converges 
to a function f at each point of D. Then f G C a (D), \\f\\ < M. 

Proof. Since \f m \ + (2R) a H a (f m ) < \\f m \\ < M, we also have 

\f m (x) - f m {x')\ < (2R)- a M\x - x'\ a 
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for all x,x' G D. Letting m — > oo in the above inequality, we conclude 

\f(x)-f(x')\<(2R)- a M\x-xT, 

and 



Therefore II f II < M 



|/| + {2R) a H a {f) < M. 



Lemma 2.7. Let {f m } be a sequence in Cq +c *(D), with ||/ m ||^ < M, and if {d^fm}, for all 
f3, \/3\ = k, are Cauchy sequences in the norm | • • • |, then there is a function f G Cq +01 (D) such 
that \d p f m - dPfl ^ asm ^ oo for all f3, < \/3\ < k, and with \\f\\^ < M. 

Proof. For (3, \/3\ = k, consider = d 13 f m , then G C a (D) and \\g!^\\ < M. Applying to 
Lemma 2.6, we have functions g 13 G C a (D) such that \g^ — g^\ — > as m — > oo, with Wg^W < M. 
Define / by 

f(x) = f r 1 ■■■ r \y. g^^Xdtdh . . ■ 

Jo Jo Jo V m=k J 

Since / G C* +a (D), we have 

f m (x) = \ i f ki ■■■ r & f m (tx)dtdt x ■ ■ ■ dt k A x ? 

\rn-i. I Jo Jo Jo J 



Therefore 

•1 rt k 



fm(x) - f(x) = J2 { l l fc " ' • • l 1 \&fm{tx) - g^tx^dtdh ■ ■ • dtfc-xja^, 



whence 



l/m-Zl^E^-A 



k\ 

\B\=k 



which goes to as m — > oo, implying f m — > f in the norm | • • • |. 

For I7I = / < k — 1 we want to prove that {d 1 f m (x)} are Cauchy sequences. Indeed, Since 
f m vanishes up to k — 1 order at the origin, then d 1 f m vanishes to k — 1 — / order at the origin. 
Thus, we have the formula, 

/l rtk-i-i fjk-l 

J dlM^ 'f m ( tx ^ dt ' ' ' dt k-i-i 



J2 & ) 8nf m {tx)jPdt---dt k - 1 - l 

\p\=k-i 

v r ••• r~ i_i dPerf m {tx)dt---dt k . 1 . l x^ (11) 

,-, =fc _,^o Jo 
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Then 

dV m {x) - 8T>f m .{x) 



= f ■■■ T' 1 '' d^{f m (tx)-f m/ (tx)}dt---dt k _ 1 _ l x^, (12) 

\p\=k-i Jo Jo 

Then 

\&f m (x) - dy m ,{x')\ < \^fm ~ Pfmfl (13) 

[ >' \P\=k 

which proves that {d^ f m {x)} (\^\ < k — 1) are Cauchy sequences since {d& f m (x)}(\/3\ = k) are. 
We assume that for = I < k — 1, d 13 f m (x) converges to g' 3 in norm | • • • |. Thus, applying 
Lemma 2.2, we have 

U^-j^l^d'Ux^x'-xf 

1=0 ' \8\=l 



k\ 

< ^-{2Ry a \\f m \\^\x' - x\ k+a < — (2R)- a M\x' - x\ k+a , 
k\ k\ 

which is independent of m. Letting m — >■ 00, we have 

W) - m - J2 \ E 9^ x ' - x ^ ^ 2 ^m- a M\x' - x \ k+ «, 
1=1 ■ w\=i 

which implies, by definition of differentiability, = d@f. This implies ||/||^ < M by taking 
limit from ||/ m ||^ < M. The convergence for \/3\ < k — 1, follows from that of = k by 
(12). I 

Lemma 2.8. The function space Cq +ol (D) equipped with the function || • • • ||^ is a Banach 
space. 

Proof. Let {f m } be a Cauchy sequence in || • • • \\( k \ Then for any e > 0, there is m so that if 
m, m! > m it holds 

||/ m -/ m '|l (fc) <e, (14) 

which implies 



. 11 (*) _ II f ,||(*)| < e 

/mil 1 1 J m' 1 1 I ^ c i 

which implies {||/m||^} are Cauchy sequence and therefore bounded by say M. Also by defi- 
nition of || • ■ • \\( h \ (14) implies d 13 f m for \(3\ = k are Cauchy sequences in | • • • |. By Lemma 
2.7, there is a function / e C fe+Q (D) such that \d^f m - d^f \ ->■ for |/3| = fc as m ->■ 00, and 
(*) < M. 
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Now the sequence {f m — /m'}m'=m+i ^ bounded in || • • ■ ||^) by e, and converges to f m — /, 
with 

d P Um-fm')^d\f m -f) 

in | • • • | as m' — > oo, \(3\ = k. By Lemma 2.7 again, it holds 

\\fm-f\\ {k) <e, 

which implies that f m — > f in || • • • ||^. The proof is complete. I 

3 Newtonian potentials and Holder estimate 
3.1 Definitions and basic facts 

The fundamental solution of Laplace's equation is given by 

r(a-y)= , 9 1 s \x-y\ 2 ' n 
n{2 — n)u n 

for n > 2, 

T(x -y) = — ln|x - y\ 

for n = 2. For an integrable function / on D, the Newtonian potential of / is the function 
N(f)(x) defined on W l by 

W)(x)= / T(x-y)f(y)dy. 



J D 



In this paper, our point view is to consider A/" as an operator acting on function space. The 
following is well-known. 

Lemma 3.1. Let f e C a (D). Then it holds, for x G M(D), 

daN{f){x) = [ dtjFix - y)(f(y) - f(x))dy - ^f(x) 

J D n 

The proof of this lemma can be found in [F, p. 204] (also also [GT]) and is based on the 
following (not so obvious) integrals: 

for n > 3, x G Int(-D) and 

m)^) = -\\A 2 + \R 2 i}o g ^+\) 

for n = 2, x G Int(D). 

Now we discuss a technical result in order to deal with Holder estimate of functions for 
operator M . Let x be an interior point of D, and D the intersection of D with the open ball of 
radius p and center x. Namely 

D = Dn{y eR n \\y-x\ < p}. 

The following is essentially proved in [F] and fundamentally important in our approach that 
follows. 
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Lemma 3.2. There is a constant C only depending on n (independent of R and p) such that 
for 1 < i, j < n 



j d i:j T(x - y)dy 

J D\Dq 



<C,xe Int(D). 



Remark 3.3. In fact, Lemma A. 14 in [F, p. 207] assumes that p < j. But the proof would go 
through essentially without the restriction. The reason is that Remark 4 in [F, p. 209] suggests 
in order to estimate Holder constant it suffices to consider the case \x — x'\ < j say. Otherwise 
when \x — x'\ > j, one would have 

\f(x)-f(x')\ < 2supl/(x)l 



:f) Q 



But this estimate would be very bad for our consideration in which we would take R — > 0. 

We note Lemma 4.4 in [GT] states that M maps C a (\D) to C a {D) continuously. In fact, 
more is true. The following result is essential to our method and seems to have been overlooked; 
and we give here a complete proof using lemmas above. 

Theorem 3.4. If f e C a (D), then M(f) G C 2+a (D). Moreover, there is a constant C(n,a), 
independent of R, such that 

\\Af(f)\\ i2) <C(n,a) 



Proof. First let us recall the norm 

\Mf)\\ l2) = max {||^(/)||}. 

l<i,j<n 

Let 

<j>{x)= [ d ij r(x-y)(f(y)-f(x))dy. 
Then, d i:j N{f) = <f>(x) - S -^f(x), and therefore 

\\Wf)\\<U\\ + 
Now we proceed to estimate ||0||. First, for x G Int(i^), 



(x)\ < [ \d ij Y{x-y)\\f{y)-f{x)\dy 

JD 

Jd \x-y\ n 



< CH a [f] / -^drdS = C(n,a)H a [f]R a - (15) 

JS"- 1 Jo r 

To compute the Holder constant of <j>, let x,x' be two (distinct) points of D. Let B(x,p) be 
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open ball of radius p = 2\x — x'\ with center at x. Let D = D fl B(x, p). We consider 



(j)(x) — 4>(x') 



D 

r 

D\D. 

r 

D\D, 



d ij T(x-y)(f(y)-f(x))dy- / dM - y)(f(y) - f(x'))dy 



D 



d l3 T(x - y)(f(y) - f{x))dy + f d l3 T(x - y)(f(y) - f(x))dy 

d ij T(x'-y)(f(y)-f(x'))dy- [ d i3 T{x' - y)(f(y) - f(x'))dy 
o Jd 

= [ (d ij r(x-y)-d ij r(x'-y))(f(y)-f(x))dy 
Jd\d 

+ (f(x')-f(x)) [ d l3 T(x'-y)dy 
Jd\d 

+ / d i3 T{x-y)(f(y)-f(x))dy- I d i3 T{x' - y)(f(y) - f(x'))dy 



D 



Do 



= h+h + h + h. 

We are ready to estimate each of I 3 : 



(16) 



\h\ 



D\D 



(d i3 r(x -y)- d i3 V{x' - y)){f{y) - f(x))dy 



< H a [f]\x-x'\ [ \Vd l3 T(x-y)\\x-x'\ a dy 

Jd\d 

< CH a [f]\x-x'\ [ \x-x'\ a \x-y\- n ~ l dy 

JD\Do 



'D\D 
»2R 



r a r -n-l r n-l 

= CH a [f]p a . 



dr 



(17) 

where x is a point on the line segment between x, x', and a polar coordinate is used at x, for 
which we have for y G D \ D , p < \y — x\ < 2R. 

By Lemma 3.2, which is a crucial difference from [GT. Lemma 4.4], we have 



\h\ = 



(f(x')-f(x)) [ d l3 T(x>-y)dy 

JD\Do 



< CH a [f]\x-x 



l\a 



(18) 



Similarly, 



/ d i3 r(x-y)(f(y)-f(x))dy 

J Do 

< H a [f] [ \d t3 T(x-y)\\x-y\ a dy 

< CH a [f] [ \x-y\- n \\x-y\ a dy 

J Do 



= CH a [f}\x-x'\ a - 



(19) 
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For I4, the estimate is identical to I3. Combining (15)-(19), we complete the proof. I 

4 An integral system 

We consider the integral system of (4) 

u = h + J\f(a). (20) 

Namely, 

Ui = hi+ M(ai) 
u 2 = h 2 +Af(a 2 ) 

Un = hN + N((In) 

where h = (hi(x), hjy(x)) with hj(x) being any harmonic function, and Af(a) is the Newtonian 
potential of a, namely 

Af(a i )(x) = [ T(x-y)a\y,u(y),Vu(y),V 2 u(y))dy. 

It is clear that any solution of (20) is a solution to (4). We further modify the equation (20) to 
fit our Banach space C* +a (D). To do so, we let, for any / G C 2+a (D), i = 1, N, 

u*(f)(x)=M(a\y,f(y),Vf(y),V 2 f(y))(x) (21) 

and define 

4 (/)(x)=^(/)(x)-^(/)(O)-^^(^(/))(O)^-- d^WMx^. (22) 

j=l k=il;k,l=l 

First we remark that the last term subtracted is harmonic since k ^ I. Therefore AO*(/) = 
Au/(/). We note that by Theorem 3.4, if / e C 2+a (D), then u\f) G C 2+a (D), and e*(/) G 
Cq +oi (D) by construction (22). More importantly, we have 

^(6 i (/))(0) = (23) 

if k ^ I. Now we introduce the Banach space from which we will seek our solutions. Define 

B(R) = C 2+a (D) x • • • x C 2+a (D) 

which consists of A^-copies of Cq +01 (D), and define the function on B(i?) as 

By Lemma 2.8, (B(i2), || • • • H'- 2 ^) is a Banach space. Now we are ready to consider a map as 
follows 

: B(R) -)• B(R) 

e(f) = (e\f),...,e N (f)). 

In order to apply a fixed point theorem on the Banach space B(i2), which has R as a parameter, 
we will have to estimate ||6(/) - Q(g)\\ (2) , and ||6(/)||( 2 ) for any f,g G B{R). This is to be 
done in the next subsections. 
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4.1 || 9(/) - Q(g) || (2) estimates 

It suffices to estimate ||0*(/) - 0*(^)|| C2) - To this end, we see from (22) 

n 

\\®V) - 0G?)ll (2) < - ^G?)ll (2) + E l<W(/) - «^))(o)|. 

k,l=l 



(24) 



First we have ||u/(/) - ^G?)ll (2) = \\Af(a l (; f, V/, V 2 /) - a l (-,g, Vg, V 2 g)) \\& by (21). By 
Theorem 3.4, we have \\J\f(4>) \\^ < C\\4>\\ for any <fi E C a (D), where C is a constant dependent 
only on n, in particular, independent of the radius R. Therefore we have 

\\co\f) - co\g)\\W < C|K(-, /, V/, V 2 /) - a\-, g, Vg, V 2 g))\\ . 

Now we are estimating ||a l (-, /, ...) — a l (-,g, ...)||. First we use coordinates for x = (xk),p = 
(Pj)iQ = Wk)> an d r = (rL)- F rom this point on, we will use constant C depending only on 
n, N, a, which varies from line to line. We begin with 



a\x,f,Vf,V 2 f)-a\x,g,Vg,V 2 g) 

^-a\x, tf + (l- t)g, tVf + (1 - t) Vg, tV 2 f + (1 - t) V 2 g)dt 

N n N N n 

E Mf: - gj) EE &Mfi - 9j) + E E ciidkWi - 9j) 

3=1 k=l 3=1 j=l k,l=l 



where 



a\x, tf + (l- t)g, iV/ + (1 - t)Vg, tV 2 f + (1 - t) V 2 g)dt 



1 _d_ 

^-a\x, tf + (l- t)g, tVf + (1 - t)Vg, tV 2 f + (1 - t)V 2 g)dt 
o dq{ 



1 d 



o dr 3 



a l (x, tf + (l- t)g, tVf + (1 - t)Vg, tV 2 / + (1 - t) V 2 g)dt. 



ki 



Taking norm || • || on (25), and using algebraic property of the norm, we obtain 



N n 



(25) 



(26) 
(27) 
(28) 



\\a\x,f,Vf,V 2 f)-a\x,g,Vg,V 2 g)\\ 

N n N 

< E W A Mf 3 -9i\\ +EEll^HIIW- -9i)\\ +E E \\<%i\\\\*Wi -9:)\\ 

3=1 k=l j=l 3=1 k,l=l 

N n N N n 

< ii/-5iiEii^n + ii/-^ (1) EEii^ii + ii/-^i (2) EEii^ii 

3=1 k=l 3=1 3=1 k,l=l 

N n N N n 

< C(R 2 EKII+^EE KW + E E II^IDII/ - 9\\ (2) (29) 

3=1 k=l 3=1 3 = 1 k,l=l 
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where we have used, according to Lemma 2.5, that 

||/ - g\\ < CR 2 \\f - s|| (2) , and ||/ - g\\W < CR\\f - g\f\ 

Since we will apply Fixed point theorem for a closed subset of B(i?), we consider the following 
closed subset 

A(i?, 7 ) = {/GB( J R)|||/||( 2 )< 7 }. 

Let 

W k = tV k f{x) + (l-t)V k g{x), 
for k — 0, 1, 2. We need to study the range of W k for /, g e A(R, 7). The following is needed. 
Lemma 4.1. If f,g e A(i2, 7), £/ien 

|W fc | < Ci? 2 " fc 7 

/or fc = 0,1, 2. 

Proof. This follows from Lemma 2.4 and 2.5 since |W fc | < ||V fc /|| + ||V fe (?||. H 

We now consider the compact set 

E(R,j) = Dx{pe R N \\p\ < Ci? 2 7 }x{g G 1R ^ 1R™ 1 1 g | < CR^}x{r e Sym(n)®M n ||r| < C 7 } 

In order to carry out estimates of (29), we need to introduce constants that would measure 
solvability of partial differential equations we consider in this paper. 

A[R,j[ 

H*[R,l] 
B[R, 7] 
H*[R,l] 

c\Rn] 
h c ;\R,i] 

Now we need more important constants on Lipschitz in variables of r. Here we denote the 
Lipschitz constant in r as follows 



= max{ 



da 1 



d Pj 
max{H a 
da 



\i = l,...,N;j = l,...,n} 



E(Rr/) 
da 1 



d Pj 



\i = 1, AT; j = 1, ...,n} 



= max{ 



dqi 



E(R,j) 

\i = l,...,N;k,j = l,...,n} 



ma.x{H a 
da 



E(R,-y) 

da*' 



\i = l,...,N;k,j = l,...,n} 



= max{ 



dr li 



\i = 1, ...,N;j = 1, ...,n} 



= m&x{H a 



E(R, 7 ) 

da 1 



dr li 



\i = l,...,N;j = l,...,n} 



E(R,j) 



(30) 
(31) 
(32) 
(33) 
(34) 
(35) 
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Now we define: 



Ht[B,j\ 






= max{i/x 


H?[R,i] 


= max{Hi 



%JB(i?, 7 ) 

da 



da 1 
dr li 



i = l,...,N;j = l,...,n} 
\i = l,...,N;k,j = l,...,n} 
\i = 1, ...,N;j = 1, ...,n} 



4.1.1 Estimate of || 

It is obvious that 

|^-| < A[R n }. 
Now we are estimating H£[Aj\. To do so, we begin with 

C 1 (9(7* 

AjW-Ajix') = / ^(x,W°(x),W 1 (x),W 2 (x))dt 
Jo VPj 

f 1 da i 

Jo dPj 
da 1 
d Pj 



-f 

Jo 

f 1 da i 
Jo dPj 



+ 



+ 



1 da i 

o dPj 
1 da 1 

o dpj 
1 da^ 

o dPj 
1 da i 

o dPj 



x',W\x'),W\x'),W 2 {x'))dt 
x,W°{x),W 1 {x),W 2 {x))dt 
x',W°(x),W 1 (x),W 2 (x))dt 
x',W°(x'),W 1 (x),W 2 (x))dt 
x',W°(x'),W 1 (x'),W 2 (x))dt 
x',W°(x'),W 1 (x'),W 2 (x))dt 
x\W°{x'),W\x'),W 2 {x'))dt 



and it follows 



\A 3 (x) - A 3 (x')\ 

< H£[R, >y]\x -xT + H£[R, 7 ] - f j (x')\ + \g j (x) - g> 



N 



3=1 



N 



+ H^R^} £(|V/'" (*) - Df(x')\ + \Vg j (x) - Dg\x f )\Y 



N 



+ H£[B,j\ ]T(|V 2 f (x) - V 2 /V)l + |vV(s) - VV(^)I)- 
j'=i 
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We remark here that C 2 regularity is used to Lipschitz estimate of r variables, and if a is 
independent of r, the regularity of C l+a is enough for the estimate. This fact will be used in 
proving Theorem 1.2 and Corollary 1.4. 

Now we need a lemma on Lipschitz property of C™ +a (D). 

Lemma 4.2. Let f e Cq +01 (D). For x,x' G D, and \P\<k- 1, we have 

\&f{x) - dPf(x')\ < ||/|| (l,3|+1) |x - x'\. 

Proof. We have 

d^f(x)-d^f(x') = £ j t d?f(tx + (l-t)x')dt 

= / Vd p f{tx + (1 - t)x') ■ (x — x')dt. (41) 
Jo 

It follows from (41) that 

\dPf(x)-dPf(x')\ < \\f\\m^\ x - x '\. 

I 

Applying Lemma 4.2 and 2.5, we have 

\f(x) - f(x>)\ < \\f\\W\x - x'\ < 3nR\\f\\^\x - x'\, 
\d j f{x)-d j f{x')\<\\f\\^\x-x'\. 

It follows from (40) 

\Mx)-Mx')\ 

< H*[R, 7 ]|x -x'\ a + H£[R,>y]N(3nRr(\\f\\<n + \\g\\ {2) T\x - x'\ a 
+ H^[R n }N(\\f\\^ + \\ g nx-xT 

N 

+ Hf[R n ]J2(Ha[V 2 f j } + Ha[V 2 g j ])\x-x'\ a . (42) 
j'=i 

Now we have by (42) 

\A j (x)-A j (x')\ < (#f[#,7]+2^(3ni?) a 7 Q /^[^ 
which implies 

H a [Aj\ < H^[R, 1 ]+2N(3nRr 1 a H^R n ]+2N 1 H^[R n ]+2N 1 Hf[R n ](2R)- a (43) 

and 

\\AA\ = \A J \ + (2R) a H a [A 1 ] 

< A[R, 7] + (2R) a (l + 2N(3nR) a ~f a + 2Nj)H£[R, 7] + 2N 1 H^[R, 7]. (44) 
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Similarly, we have estimates 

H^ll = \B 3 \ + {2RTH a [B{} 

< B[R, 1 ] + (2R) a (l + 2N(3nR) a 1 a + 2N 1 )H^[R, 1 } + 2N 1 H^[R, 1 ], (45) 

and 

\\CL\\ = \C{ l \ + {2RTH a [C{ l ] 

< C[R,j] + (2R) a (l + 2N{3nR) a 1 a + 2N*y)H%[R,i\ + 2N>yH?[R,i\, (46) 

To simplify the notation, we denote the right side of (44), (45), and (46) respectively by 
8 A (R,^),8 B (R,^) and 8 C (R,^). Then by (29) we, have 

\\a\x,f,Vf,V 2 f)-a\x,g,Vg,V 2 g)\\ 
< C(n, N, a)(R 2 8 A (R, 7 ) + R8 B (R, 7 ) + 5 C (R, 7)) 11/ - <?|| (2) - (47) 

Lemma 4.3. Let f e C a (D). Then it holds 

l^(/)(0)|<C(n,a)||/||. 

Proof. By Theorem 3.4, we have 

IW)(0)l < l^(/)l < IIW)|| (2) < C(n,a)\\f\\. 

I 

Finally, we combine (24), and Lemma 4.3 to conclude that 

\\Q\f)-&(g)\\ (2) <S(R,i)\\f-9\\ m (48) 

where 

5(R }1 ) = C(n,N,a)(R 2 5 A (R, 1 ) + R5 B (R, 1 ) + 5 c (R,i)) (49) 

5 A (Rn) = A[R, 1 ] + (2R) a (l + 2N(3nR) a 1 a + 2N 1 )H^[R, 1 ]+2N 1 Hf[R, 1 ] (50) 

8 B {R,l) = B[R n ] + (2R) a (l + 2N(3nR) a 1 a + 2N 1 )H^[R, 1 }+2N 1 H^[R, 1 ] (51) 

Sc(Rn) = C[R, 1 ] + (2R) a (l + 2N(3nR) a 1 a + 2N 1 )H^\R, 1 ] + 2N 1 H^[R, 1 }. (52) 

4.2 ||e(/)||( 2 ) estimates 

Our proof here is mostly similar to that of (48). We will point out minor differences without 
repeating the arguments. It suffices to estimate ||O l (/)||^- To this end, we see from (22) 

n 

\\&(f)\\ {2) < lk(/)ll (2) + E M W))(0)|- (53) 

k,l=l 

First we have || < 2 > = ||./vV)||( 2 ). By Theorem 3.4, we have ||WV)II (2) < C||a*||, where C 

is a constant dependent on only on n, in particular, independent of the radius R. 
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A 3 = 


9 i, 
Jo dPj 


[x,tf,tVf,tV 2 f)dt 


B{ = 


Jo dqi 


[x,tf,tVf,tV 2 f)dt 


C L = 


r a .o- 

Jo dr 3 kl 


(x,tf,tVf,tV 2 f)dt 


Dj = 


C 9 a*. 
Jo 9xj 


[x,tf,tVf,tV 2 f)dt. 



Now we are estimating ||a l ||. First we use coordinates for x = (xj),p = (pj),q = (ql), and 
r = (r J kl ). We begin with 

a\xJ,Vf(x),V 2 f(x))-a l (0,...,0) 
^-a\tx,tf,tVf(x),tV 2 f(x))dt 

o 

n N n N N n 

= E D i x i + E mm + E E &Mfi) + E E <W(/;) ( 54 ) 
j=i j=i fc=i j=i j=i fc,/=i 

where 

/•i « 

(55) 
(56) 
(57) 
(58) 

We notice here we have kept notations similar to (26)-(28), although the integrand is different. 
By virtue of the same argument as for (24), we can arrive to the following estimate 

\\a\x,f,Vf,V 2 f)\\ < 1^(0)1 + 3^11^11 

N n N N n 

+ C(i? 2 ^K.|| +J R^^||^||+^^||^||)||/||( 2 ). (59) 

j=\ k=l j=l j=l k,l=l 

Similarly, we need to define constants D[R, 7)], H®[R, 7], and H®[R, 7] as in (30)-(38). With- 
out repeating, we will come out estimate 

l|e*(/)|| (2) < v(R,i) (60) 

where 

•n(R,i) = C(n, N, a)(\a(0)\ + R5d{R, 7) + 7)) (61) 

6 D (R,*y) = D[R, 1 } + (2R) a (l + 2N(3nR) a 1 a + 2N 1 )H°[R, 1 } + 2N 1 H?iR n ], (62) 

and 8(R, 7) is given by (49). 
4.3 A general estimate 

Here we collect the estimates together for a later quick reference. 

Theorem 4.4. Let : B(R) ->■ B(R) be defined as in (22). Iff,ge A{R,~f), then 

\\e(f) - e(g)\\W < 5(R, 1 )\\f-g\\^ (63) 
l|e(/)|| (2) < r/(i?, 7 ) (64) 
where 5(R,j) and n(R,j) are defined by (49), (61) respectively. 
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5 Proof of theorems 



In this section, we will give proofs of all results presented in the introduction. It suffices to work 
with a of class C 2 or C 1+a since the regularity of A insures higher regularity. 



5.1 Non-radial functions 

In order to construct solutions that are not radial, it is helpful to have the following simple 
lemma for this purpose. 

Lemma 5.1. Let u be any function of C 2 (D). If V 2 m(0) is not XI nX n for some A G R, then u 
is not radial. 

Proof. If u is radial, then there is a function v such that u(x) = v(r) where x = ruj,u G S" 1-1 . 
Obviously v is C 2 for r > 0. In fact, we have, for r > 0, 

k,l=l 

Taking limit on both side of the equation above, we have 

n d 2 u 
]imv"(r) = Y^lim- — —(x)u k ui 

r-»-0 ^ x^O dXhOXi 

k,l=l K 1 

n o2 

^ dx k dxi UJkL ° 1 ' 
k,l=l ft 1 

In particular \im r ^ Q v"(r) exists and we denote it by A, and V 2 w(0) by A. We conclude that 
A = uiAui 1 - for any u> G namely, a; (A — A/)^ = for any u G S n_1 . This means that 

A = XI, a contradiction. S 



5.2 Proof of Theorem 1.1 

Here we will prove a result slightly more general than Theorem 1.1. 

Theorem 5.2. Let a(x,p, q, r) = (ai(x,p, q, r), a^{x,p, q, r)) be of class Cf oc (2 < k < oo, < 
a < 1), where x G W l ,p G R N ,q G IR™®]^^, and r G Sym(n) ® R^. There is a (small) constant 
5 depending on n, N, a such that if 

a(0) = 0, (65) 
|V r o(0)| + |V 2 a(0)| < 5, (66) 

then the following system: u(x) = (ui(x), ...,un(x)) : {\x\ < R} — > M N , 

Au(x) = a(x,ti(3;),Vii(x),V 2 n(x)) (67) 

has infinitely many solutions of C k+2+a (D) of vanishing order two at the origin for sufficiently 
small values of R. Furthermore these solutions are not radially symmetric. 
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In the proof, our goal is to find R, 7 sufficiently small so that we have 

v(R,i) < \ 

Once we have these, we can consider an operator T defined by 

T(f) = h + e(f) 

where / G A(R, 7) and h is a harmonic vector-function such that h(0) = V/i(0) = 0, and 
WhW^ = 2. Therefore T is a contration operator from A(R, 7) to A(i2, 7), for which there is a 
fixed point that becomes a solution of (20). 

To this end, we first simply S(R, 7), r}(R, 7), and we can replace them by 

S(R n ) = C(C[ J R,7]+7^f 7 [^7])+^,7) (68) 
17(^,7) = C(|a(0)|+7C[ J R,7]+7 2 ^f[ J R,7])+^,7)- (69) 

where C is a constant only depending on n,N,a, and 7) is such that lim^o 7) = 
for each 7 > 0. We now give estimates of C[R, 7], iff 7] in terms of conditions (l)-(3). Let 
a - = be a component variable of r. We want to estimate the Lipschitz constant of d a a. In 
fact, we have 

d a a l (x,p,q,r) - d a a l (x,p, q,r') 
1 d ■ 

—a l (x,p,q,tr + (l — t)r)dt 
dt 

= /E^^-)^-^)*- (70) 
Here a l (-) is defined naturally as given in the integrand. It follows from (70) that 

H\[d a a % ]\ E (R, 1 ) < \V 2 r a l \ E (R rt ) 

= \V 2 r a i (0)\ + o(R + 1 ), (71) 

where o(R+ r y) — > as R, 7 — > by continuity of C 2 smoothness of a and lim^ ^o E(R, 7) = {0}. 
From (71) and definition (38), we have 

H^[R n ]<\V 2 r a(0)\ + o(R + 1 ). 

On the other hand, we have 

d a a l (x,p,q,r) — d a a l (0) 

f 1 d ■ 

= J —a\tx,tp,tq,tr)dt 

+ ( 72 ) 
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Notice that for (x,p, q, r) G E(R, 7), we have \x\ < .R, |p| < CR 2 ~f, \q\ < CR'y, and |r| < C7. 
Hence by (72), putting terms with R,j together, we have 

\d„a l \ E{Rn) < C(|V r a*(0)| + | V^(0)|) + e(R, 7) 

where 7) — >• as R — > for each given 7. This implies 

C[R,j] < C(|V r a(0)| + |V r 2 a(0)|) +e(R,>y) 

Letting 

r = |V r a(0)| + |Vj?a(0)|, 

we can have 

<J(i2, 7) = C(t + 7(t + o(i2 + 7)) + 7) 
7/(i2, 7) = C( 7 r) + 7 2 (r + o(i? + 7)) + e(R, 7). 
By o(R + 7) — >■ 0, there exist i? , 7o(< 1) such that 

o(i2 + 7) < Ct 

for < Rq, 7 < 7 . Hence we have 

W To) <Ct + 7o 2Ct + e(i2, 7o ) 

v(R,lo) < 7o3Ct + £ ( j R, 7o ) . 
Now choose <5 = If r < <5, we have 

S(R,lo)<^ + ^ + e(Rno) 

v(Rno)<^ + e(Rno)- 

Finally, we choose R small so that e(R, 70) < min {^, ^}. It follows that 

S(R,io)<^ 

This is equivalent to that for /, g e A(i?, 70) 

l|0(/)-0(^)ll (2) < ^ll/-^ll (2) (73) 
ll©(/)ll (2) < f- (74) 

Now we are ready to apply Fixed point theorem for R, 70 . Let h be any harmonic homogenous 
polynomial of degree 2 so that ||/i||^ 2 ' ) = y. Here let us be more specific. Let h = Ylki=i a ki x kXi 
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where a k i € K , = a«fe. It is easy to see that /i is harmonic if and only if the trace of a is 
zero, i.e., Ylk=i akk = i®}- We a ^ so see ^ na ^ 

||/i|| (2) = max \a kl \. 

l<k,l<n 

So we take a harmonic polynomial of degree 2 such that < max.i< k ,i<n \o>ki\ < -f- Then we 
consider the map 

T(u) = h + e(u) 

which maps A(R, 70) to A(i2, 70) as contraction map by (73), (74). So T has a fixed point u in 
A(i2, 70). We claim the vanishing order of u at the origin is 2. In fact, by (23), 

d k d lU (0) = d k dM0) = a kl ^ {0} 

for k 7^ I, and also w(0) = and d k u(0) = by the construction. At the same time, for different 
{a k i} we have different solutions u. Now we prove all these solutions are not radial. Assume 
there is i such that (a k °) 7^ {0}. If (a]j.°) = XI for some A G K, then since the trace of (a^°) is 
zero, we conclude that A = 0, which implies (a^°) = 0, a contradiction. By Lemma 5.1, u l ° is 
not radial. So u is not radial neither. This completes the proof. 

Remark 5.3. Given 70 as in the proof above, we actually prove that the solution space can be 
parameterized by at least tt-^N parameters from the coefficients of the harmonic polynomials 
of degree 2. In fact, the different choice of a ki produces different u which is defined in the same 
domain of radius R and with different Hessian at the origin since d k diu(0) = d k dih(0) = a k i by 
(23). Therefore there exist infinitely many solutions of vanishing order two at the origin. This 
remark also applies to other theorems 

5.3 Proof of Theorem 1.2 

First we consider the case where c = 0,(7 = as in (6), (7). Since a is independent of r, we 
have C[R, 7] = H^[R, 7] = H^[R, 7] = 0. Here we note that C 1,a regularity of a is only needed 
with a careful inspection of the proof of Theorem 4.4. Therefore we can replace constants as 
follows 

S(R, 1 )=e(R, 1 ) 
77(i?,7) = C|a(0)|+e(i2, 7 ) 

where 

lim e(R, 7) = lim e(R, 7) = 

for each 7 > 0. Now we choose 70 large enough that ^ > C|a(0)|. Then we choose R sufficiently 
small that e(R, 70) < \ and e(i?, 70) < -j. As a result, we have 

S(R,io)<\ 

As in Theorem 1.1, we find a solution in A(R, 70) which vanishes up to order 1 at the origin. 
To get general case, we consider a new system 

Au = a(x, u + Co + ci • x, V(-u + c + c\ • x)) = a(x, it, V-u) 
We can solve this system for u. Then u = u + cq + c\ • x is the solution we are seeking for. 
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5.4 Proof of Theorem 1.3 

Here we fix R, and choose 7 small to prove the existence of semi-global solutions. Since a is 
independent of x, so we have D[R, 7] = H^[R, 7] = H®[R, 7] = 0. Therefore we can replace 
5(^,7)^(^,7), using (8), (61), by 

#, 7 ) = W 7 ) (75) 
where 5(R, 7) is still given by (49). It suffices to prove that for each given R, we have 

lim 5(R, 7) = 0. (76) 

If this is proved, then we can choose 7 so that 

v(R,io) <|- 

Our next goal is to show (76). Indeed, since a is independent of x, we will take E(R, 7) as 
E{R,i) = {pe R N \\p\ < CR 2 ^} x{qe R N ®R n \\q\ < Ci?7} x {r G Sym(n) <g> R n | |r| < C7}. 



We notice that as set, 



lim£(i?, 7 ) = {0}. 

7— >o 



This is important for what follows in proving (76). Let a be one of component variables of 
{p, q, r}. We have by condition (9) 

f 1 d 

d a a l (p,q,r) = J —a\tp,tq,tr)dt 

= ^ $~>A^(-)^ + Y, d <tMU& + E^A^KA ( 77 ) 

Since \pj\ < Ci? 2 7, \q J k \ < CR^y, and \r J kl \ < C7, we have from (77) 

\d a a % \E{R, 1 ) < C\\a\\ C 2(E{R,j))l, 
which implies, by definitions (30), (32), (34), 

A[R,i\,B[R,i\,C[R,l\ < C\\a\\ c , mi))1 . (78) 

Here, of course, C 2 (E(R,^)) denotes the maximum norm of V 2 a on the set E(R, r y). On the 
other hand, we have 

d a a\p,q,r)-d a a\p\q'y) 
f 1 d • 

= J —a l {tp+{l-i)p,tq+{l-t)q,tr + {l-t)r')dt 

= J*T, d pM(')(Pj-Ji) + £^6U(.)(^ - q'i) +J2 d ri l d « ai (-)( r li - r'L)dt. (79) 
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From (79), we have, by definition (36), (37), (38), 

H£[R,>y],Hj![R,>Y],H£[R,>y] < C||a|| ca(£(JI>7)) 7 1 - a . (80) 

Similarly, from (79), we have 

H?[R,i\,H*[R, 7], flf [i2,7] < C\\a\\ cHE{Rn)) . (81) 

Finally, we substitute (78), (80), and (81) into (50), (51), (52) and (49), we see that 5(R,i) is a 
function in 7 and 7 1 ~ a , which proves (76). The rest of proof is similar and we omit it. 



5.5 Proof of Corollary 1.4 

It follows from Theorem 1.3 immediately since C 1+a is only needed due to the remark in the 
proof of Theorem 1.1. 
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